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IHTRODUCTION 

Through  simulation  the  performance  of  an  orcanizatlon  or  some 
part  of  an  orcanization  can  be  represented  over  a  certain  period 
of  time  (Martin,  I96I).  This  is  accomplished  by  devising  mathe- 
matical models  representing  the  components  of  the  organization 
and  'oroviding  decision  systems  for  these  components  to  represent 
the  various  interrels-tionships.  A  digital  computer  may  be  used 
to  progress  through  the  model  for  various  intervals  of  time,  paus- 
ing at  the  end  of  each  interval  to  compute  the  interactions  be- 
ti.'cen  certain  components.   In  this  way  it  is  possible  to  obtain  a 
"simulated"  history  of  the  organization  under  certain  specified 
conditions.  By  changing  the  conditions  and  repeating  the  process 
it  is  possible  to  compare  various  parameters  or  decision  systems. 
In  -this  manner  it  is  possible  to  experiment  with  changes  in  an 
organization  without  affecting  the  actual  organization  being 
studied. 

Eandon  variations  are  characteristic  of  many  organizations. 
In  many  situations,  the  randomness  is  so  essential  that  the  com- 
mon simplifying  technique  of  using  average  values  simply  "assumes 
away"  the  problem.  VJhenever  uncerts-inty  is  s.n  essential  considera- 
tion in  an  organization,  a  method  must  be  available  to  produce 
a  similar  uncertainty  in  the  simulation  model  of  the  organization. 

The  use  of  a  computer  in  such  problems  depends  to  a  great 
c:;tcnt  on  one's  ability  to  structure  a  simulation  in  a  machine 
or  cyntcr.s  language.  Kore  specifically,  it  depends  a  great  deal 


on  the  oonfomance  of  certain  prosrammins  languases  to  tho  needs  of 
simulation  problems.  For  the  programmer  the  most  important  need 
is  that  the  language  be  problem  oriented,  since  an  organization 
will  be  reduced  to  a  mathematical  model. 

Inherent  in  all  organizations  are  many  and  varied  decisions, 
i.  e.  rules  of  action.  Decisions  may  more  accurately  be  described 
as  decision  systems  since  a  single  decision  may  be  dependent  on 
many  factors.  Decision  systems,  due  to  their  complexities,  normally 
occupy  the  greatest  portion  of  the  simulation  routing.   Thus,  the 
ability  of  a  programming  language  to  implement  decision  systems  is 
important .  ■     -.    v      \     ' 

Speed  and  efficiency  are  also  primary  factors.   Speed  refers 
to  the  time  required  by  a  compiled  decision  system  to  make  the  cor- 
rect decision.  This  time  must  be  minimized  without  occupying  ex- 
cessive amounts  of  computer  memory. 

Efficiency  refers  to  the  amount  of  computer  memory  required. 
Since  the  memory  must  provide  room  for  the  compiled  mainline 
program,  and  many,  sometimes  large,  areas  for  working  and  storage, 
the  amount  occupied  by  any  one  of  these  should  be  minimized.   The 
size  of  the  mainline  program  and,  to  a  certain  extent,  the  size  of 
the  working  areas  are  variable  and  depend  upon  the  simulation  being 
performed.  However,  since  the  subroutines  do  the  bulk  of  the  work, 
it  seems  only  natural  to  make  them  as  efficient  as  possible.   This 
is  true  especially  for  the  decision  process. 

Another  important  and  sometimes  over-looked  factor,  is  flexi- 
bility.  The  process  of  trying  several  different  conditions  or 
changes  within  the  organization  is  inherent  to  simulation.   These 


chances  are  typically  within  the  decision  process  itself.   Thus, 
provisions  should  be  nade.  to  Ha]:e  these  changes  in  the  fora  of 
data  on  the  object  level.   This  should  reduce  the  need  to  recom- 
pile a  program  each  time  a  change  must  be  made. 

Three  important  factors  in  the  selection  of  a  procraramins 
lanEuage  for  use  with  simulation  pi-oblens  have  been  discussed. 
They  are  speed,  efficiency  in  the  size  of  compiled  progrsjns  and 
their  subroutines,  and  flexibility  on  the  object  level  of  the  com- 
piled prosrams.  Another  requirement  in  Konte  Carlo  simulations  is 
the  ability  to  generate  and  use  of  random  numbers   (Hull,  19^2). 

Chance  variation  is  generally  introduced  into  systems  models 
by  probability  distributions.   The  distribution  function  may  be  a 
standard  one  (i.e.  normal,  exponential,  Poisson,  etc.)  or  an  arbit- 
rary observed  distribution  derived  from  historical  data.  The  most 
important  of  these  distributions  is  the  rectajigular  or  uniform  dis- 
tribution because  any  distribution  may  be  represented  by  random 
numbers  dratm  from  a  rectangular  distribution  by  use  of  a  probabil- 
ity integral  transformation  (Freund,  I963) . 

IBM  1620  Fortran  II  is  not  a  programming  language  that  satis- 
fies all  of  the  above  needs.   Thus  it  is  the  objective  of  this  thesis 
to  augment  the  Fortran  II  programming  language  so  that  the  needs 
specified  above  are  satisfied. 


■'■EaJidom  nujibers  generated  on  a  computer  are  called  "pseudo-random" 
numbers  and  only  appea.r  to  be  dravm  at  random  from  certain  probabil- 
ity distributions.   However,  they  are  expected  to  be  non-repeating 
for  a  "long"   sequence  of  numbers. 


?OETEAH  II 

Fortran  ("formula  translation")  II  is  a  language  that  closely 
resenbles  algebra  (IBM,  I962) ,   It  is  a  programming  language  de- 
signed primarily  for  scientific  and  engineering  calculations. 
Since  It  is  problem  oriented,  it  provides  engineers  with  a  method 
of  communication  that  is  more  familiar,  easier  to  learn,  and  easier 
to  use  than  actual  machine  language.   In  this  light,  Fortran  II 
meets  one  of  the  aforesaid  qualifications  for  a  programming  lan- 
guage that  is  suited  for  simulation  problems. 

In  addition,  Fortran  II  object  programs  are  generally  as  ef- 
ficient as  those  which  might  be  icritten  by  an  experienced  program- 
mer (IBM,  1962).   A  baiilc  of  efficient  subroutines  is  included  in 
every  Fortran  program.   These  subroutines  are  used  throughout  the 
entire  program.   The  mainline  program  consists  of  branches  to  the 
subroutines.   Since  the  subroutines  occupy  such  an  Important  posi- 
tion and  make  up  the  majority  of  the  object  program,  great  pains 
have  been  talcen  to  make  them  as  short  and  as  fast  as  possible. 

In  addition  to  speed  and  brevity,  Fortran  II  is  equipped  with 
floating  point  capabilities.  Although  this  is  not  necessary,  it  is' 
convenient  and  makes  programming  considerably  easier. 

Perhaps  the  greatest  disadvantage  associated  with  Fortran  II, 
is  the  complexity  of  the  decision  process  for  multiple  decisions. 
The  basic  element  of  the  decision  process  is  the  "IP"  statement 
which  permits  the  programmer  to  change  the  sequence  of  statement 
execution,  depending  upon  the  value  of  an  arithmetic  expression. 


The  general  form  is:  •   - 

IF(a)Hl,W2,N3 

where  (a)  is  an  expression,  and  Nl,  N2,  and  N3  are  statement  num- 
bers.  Control  is  transferred  to  statement  number  Nl,  N2,  or  K3  de- 
pending on  whether  the  value  of  (a)  is  less  than,  equal  to,  or 
greater  than  zero  respectively  (IBM,  I962).   The  "IF"  statement  will 
test  only  one  decision  parameter  at  a  time. 

As  long  as  only  one  parameter  need  be  tested  at  a  time,  the 
"I?"  statement  is  quite  sufficient.   However,  the  number  of  "I?" 
statements  required  increases  e:cponentially  with  the  number  of  para- 
meters in  the  decision  system.  As  an  example,  suppose  A  and  B  need 
to  be  tested  independently.  A  flow  diagram  of  the  process  is  as  in 
figure  1. 

This  means  that  there  are  nine  possible  configurations  of  the 
tvro  variables.   In  tabular  form  these  are  as  follows. 

A+  +  -1- 000 

B+-0  +  -0  +  -0 

Each  possible  configuration  will  require  a  seperate  branch  to 
a  seperate  or  identical  position  within  the  program.  A  program  fol- 
lowing the  logical  flow  listed  above  would  contain  the  following 
"IP"  statements: 

IF  (A)  11,12,13 

11  IF  (B)   1,  2,  3 

12  IF  (3)  ^i,    5,    6 

13  IF  (B)  7,  8,  9 


It  is  possible  to  generalize  the  nimber  of  possible  paths  as 
follows! 

Number  of  Paths  =  3  . 

vihere  n  =  the  number  of  decision  parameters.   The  number  of  "IF" 
statements  necessary  to  program  a  decision  on  (n)  parameters  would  be 

"■1-1 
i=l 
This  can  be  seen  by  sitniming  the  following  series,  representing  a 
decision  tree. 

1    2  n-1 

1  +  3  +  3  + +  3 

Thus  if  twenty  parameters  are  necessary,  there  would  be  3,^86,78^,^01 
possible  paths  and  1,7^1-3,392,110  "IP"  statements  would  be  necessary 
to  malce  the  decision  as  to  which  path  should  be  followed.   In  addi- 
tion to  this  difficulty  in  implementing  complex  decisions,  the  mach- 
ine language  programs  resulting  from  "IP"  statements  are  relatively 
slow  and  inefficient. 

Another  disadvantage  to  Fortran  II  is  that  random  numbers  are 
not  readily  available.   Certainly,  there  are  several  models  by  which 
random  number  generators  may  be  inserted  into  Fortran  programs,  but 
these  are  relatively  time  consuming  if  they  are  programmed  in  the 
Fortran  language.   The  same  models  when  programmed  in  SPS  (Symbolic 
Programming  System)  are  much  faster  and  occupy  much  less  space. 

The  choice  of  Fortran  II  for  further  study  and  implementation 
is  quite  logical  since  it  is  "easily  implemented  to  problem  solv- 
ing applications  and  still  compiles  relatively  efficient  machine 


language  programs".   It  remains,  however,  to  Implement  a  method  by 
which  decisions  may  be  made  more  efficiently,  both  in  terms  of  the 
memory  occupied  and  in  terms  of  the  time  necessary  for  execution. 
Implementation  of  an  efficient  "random  number  generator"  will  also 
be  necessary. 

■  Fortran  II  is  capable  of  using  "functions"  or  subroutines. 
Functions  are  divided  into  three  types: 

1.  Library  functions 

a.  Non-relocatable  library  functions 

b.  Relocatable  library  functions 

2.  Arithmetic  statement  functions 

3.  Fortran  subprograms    ■,  ^  •,  ' 

a.  Subprogram  functions 

b.  Subprogram  subroutines. 

The  most  elementary  type  of  function  is  the  "non-relocatable"  li- 
brary f miction.   These,  along  with  the  relocatable  library  functions, 
occupy  the  subroutine  dec)!:  of  the  Fortran  II  compiler.  This  type 
of  subroutine  is  used  for  the  Input  and  output  fxinctions  as  well  as 
the  basic  arithmetic  operations  (such  as  add  or  multiply).   They  are 
under  the  control  of  the  processor  only  and  may  not  be  called  direct- 
ly from  a  Fortran  source  program. 

Non-relocatable  library  functions  occupy  a  fixed  area  in  core 
storage.   They  will  be  present  in  all  Fortran  programs  and  in  the 
same  location  unless  the  processor  is  changed.   These  subroutines 
may  be  changed  only  with  direct  changes  to  the  processor.   It  would 
be  better  to  rewrite  the  entire  processor  than  to  attempt  to  worl: 
at  the  level  of  the  "non-relocatable"  library  function.   The  rewards 


for  direct  work  on  the  proooGsor  would  be  about  the  same  for  a  great- 
er investment  in  time  and  effort.   Hereafter,  when  reference  is  made 
to  library  functions,  it  will  be  to  the  exclusion  of  the  non- relo- 
catable variety. 

Relocatable  library  functions  will  be  compiled  into  the  object 
program  only  if  they  are  called  from  within  the  source  program  (with 
the  exception  of  the  exponential  and  logarithm  library  functions), 
and  they  will  not  occupy  a  fixed  location.   Certain  functions,  for 
instance  the  square  root  function,  are  not  always  required  and  would 
only  occupy  space  in  the  computers  memory.   These  are  included  as 
relocatable  library  functions  so  that  core  storage  is  not  needless- 
ly taken  up. 

Library  functions  are  preprogrammed  and  exist  in  a  prepared 
deck,  referred  to  as  the  subroutine  deck  of  Fortran  II.   Instead  of 
appearing  in  the  object  program  every  time  they  are  called,  they 
appear  only  once  and  only  if  they  are  called.   They  may  be  called 
by  including  the  name  of  the  function  in  an  arithmetic  statement 
within  the  source  program.  The  name  is  follovred  by  a  single  argu- 
ment enclosed  in  parentheses  which  may  be  a  constant,  a  variable 
(subscripted  or  not  subscripted),  or  an  expression. 

An  example  of  a  statement  which  will  call  two  of  the  library 
functions.  Viz,  SINF  and  SQHTF,  is: 

Y  =  A-SIWP(B*SQRTP(C) ), 

In  this  case,  the  assembled  instructions  of  the  object  program 
will: 

1.   Branch  to  the  square  root  subroutine  to  compute  the  value 


of  the  square  root  of  C. 

2.  Multiply  the  value  of  the  square  root  of  C  by  3, 

3,  Branch  to  the  SIHF  subroutine  to  conpute  the  value  of  the 
sine  of  the  product. 

^.   Subtract  this  value  from  A. 

5.   Replace  the  present  value  of  Y  with  the  value  obtained  in 
step  4, 
Only  one  value  is  produced  by  a  library  function  and  only  one  value 
may  be  used  as  an  argument,   Fortran  II  has  the  capability  of  using 
up  to  fifty  library  functions.   The  compiler  furnished  by  IBM  has 
seven.   The  process  by  which  a  user  may  add  library  functions  will 
be  discussed  later. 

An  arithmetic  statement  function  is  defined  by  one  arithmetic 
statement  and  applies  only  to  the  program  in  which  it  appears.  The 
name  of  the  function  determines  the  mode  of  the  value  that  is  com- 
puted.  The  mode  must  be  either  fixed  point  or  floating  point  as 
determined  by  the  first  letter  in  the  name  of  the  function. 

The  arithmetic  statement  function  is  defined  as  follows: 

NAME  (AEG)  =  EXPRESSION 

where  KAIiE  is  the  name  of  the  function,  AEG  is  the  argument"  and  con- 
sists of  one  or  more  variable  names  of  non-subscripted  variables; 
and  EXPRESSION  is  an  arithmetic  expression  that  must  conform  to  the 
rules  for  forming  expressions.  Examples  of  arithmetic  statement 
functions  are  as  follows: 

FIHST(X)  =  A  *  X  +  B     ■ 
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SECONr>(X,B)  =  COS(X)  *  FIHST(B) 

Y  =  SECOND(C,D). 

This  sequence  of  functions  could  be  replaced  (losing  a  great  deal 
of  versatility  but  saving  much  space  and  computing  time)  by  the  fol- 
lowing statement. 

Y  =  COS(C)  *  (A  *  D  +  B)       ■   ■ 

The  appearance  of  the  name  of  an  arithmetic  statement  function 
in  an  arithmetic  statement  serves  to  call  the  function. 

As  many  of  the  variables  appearing  within  the  expression  of  a 
function  as  desired  may  be  stated  on  the  left  hand  side  of  the  arith- 
metic statement  function  as  arguments.   The  arguments  are  really 
only  dummy  variables  so  their  names  are  unimportant  except  to  spec- 
ify the  mode.  They  may  even  be  the  same  as  names  appearing  later 
in  the  program. 

Those  variables  which  are  not  stated  as  arguments  are  treated 
as  parameters.  Thus,  if  FIHST  is  defined  in  a  function  statement  as 

FIRST (X)  =  A  *  X  +  3, 

then  a  later  reference  to  FIRST (Y)   will   cause 

A  *  Y  +  B, 

based  on  the  current  values  of  A,  B,  and  Y,  to  be  computed. 

The  arguments  of  an  arithmetic  statement  function  may  be  ex- 
pressions and  nay  involve  subscripted  variables.   Thus,  a  reference 
to 

FIRST (Z  +  Y(I)), 
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as  a  result  of  the  previous  definition  of  FIRST,  villi  cause 

A  *  (Z  +  Y(I))  +3, 

to  be  computed  on  the  basis  of  the  current  values  of  A,  3,  Y(I), 
and  Z, 

Functions  defined  by  arithmetic  statements  are  always  compiled 
as  closed  subroutines.  A  closed  subroutine  is  used  only  by  calling 
for  it  in  the  source  program  and  will  be  used  only  in  the  program  in 
which  it  was  compiled.  This  means  that  the  machine  language  instruc- 
tions are  compiled  only  once  in  the  object  progr;3un.  Calling  is  ac- 
complished either  with  the  use  of  a  "CALL"  statement  or  by  mention- 
ing the  name  of  a  subroutine  in  the  source  program. 

Fortran  subprograms  consist  of  subroutines  that  are  not  used 
frequently  enough  to  be  library  functions,  yet  due  to  their  size 
or  complexity,  they  cannot  be  defined  in  a  single  arithmetic 
statement. 

Fortran  subprograms  are  compiled  seperately  but  the  object 
program  cannot  be  executed  by  itself.   For  execution  they  must 
be  called  and  loaded  by  a  main  program.   For  this  reason  they  are 
called  subprograms.   Since  each  subprogram  is  compiled  seperately, 
the  arithmetic  function  and  variable  names  used  in  the  main  program 
are  completely  independent  of  the  function  and  variable  names  used 
in  the  subprogram.  This  means  that  very  general  subprograms  can  be 
cres.ted  and  used  with  any  main  Fortran  II  program. 

Subprograms  may  be  divided  into  two  types;  function  subprograms 
and  subroutine  subprograms.  Four  statements  are  necessary  for  their 
definition  and  use.   They  are: 
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FUNCTIOH, 

SUBROUTINE,    T  '  '";  ;        ,'1 

CALL,  and 

EETUHN, 

and  will  be  described  later. 

Although  function  and  subroutine  subprograms  are  similar  and 
are  treated  together  in  this  thesis,  thoy  differ  in  two  signifi- 
cant respects: 

1.  Function  subprograms  can  compute  only  one  value,  whereas 
subroutine  subprograms  can  compute  many  values  and  return 
them  to  the  main  program. 
•  2.   Function  subprograms  may  be  called  by  an  expression  con- 
taining its  name  but  subroutine  subprograms  are  called  by 
the  use  of  a  CALL  statement.   This  means  that  function  sub- 
programs may  be  used  in  arithmetic  statements. 
The  FUNCTIOK  statement  is  always  the  first  in  a  function  sub- 
program and  defines  it  as  such.   The  general  form  of  a  FUNCTION 
statement  is  as  follows: 

FUNCTION  "Name"  (al,  a2, an). 

"Name"  is  the  symbolic  name  of  a  single-valued  function,  and  each 
argument  (al,  a2,  ...  ,  an)  is  a  non-subscripted  variable  name, 
i.e.  must  not  be  a  member  of  a  dimensioned  array. 

In  a  function  subprogram,  the  name  of  the  function  must  appear 
in  an  input  list  or  on  the  left-hand  side  of  an  arithmetic  state- 
ment. An   example  indicating  the  use  of  a  FUNCTION  statement  and  a 
HETUEK  statement  and  which  conforms  to  this  rule  is  as  follovjs: 
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FUNCTIOS  SUK  (A, 3) 
SUM  =  A  +  B 
BETUEN. 

The  EETUEN  statement  terminated  the  function  subprogram  and 
returns  the  value  of  the  function  to  the  calling  program.  Any 
number  of  these  statements  may  be  used  but  control  must  end  with 

one  of  them. 

The  arguments  which  follow  the  function  reference  in  the  cal- 
ling program  must  agree  in  number,  mode,  and  order  with  the  argu- 
ments listed  m  the  FTOICTION  statement  in  the  function  subprogram. 
Vmen  the  argument  in  the  reference  statement  is  an  array  name  the 
corresponding  argument  in  the  PUHCTIOH  statement  must  also  be  an 
array  name  sjid  both  must  be  dimensioned  ^jithin  their  respective 
programs.   None  of  the  dummy  variables  may  appear  in  EQUIVALENCE 
statements  in  the  function  subprogram. 

The  SUBHOUTINE  statement  is  used  in  a  subroutine  subprogram 
in  the  same  manner  a  FUNCTION  statement  is  used  in  a  function  sub- 
prograjn.   The  general  form  of  the  statement  is  as  follows: 

SUBHOUTINE  "Name"  (al,  a2 ,  an). 

"Name"  is  the  symbolic  name  of  a  subprogram,  and  each  argument, 
if  any  is  specified,  is  a  non-subscripted  variable  name.   The 
SUBROUTINE  statement  defines  the  subroutine  subprogram  and  must 
be  the  first  statement  in  the  subprogram.   The  subprogram  must 
be  a  Fortran  program  and  may  use  any  Fortran  II  statement  except 
a  FUNCTION  statement  or  another  subroutine  statement, 

A  calling  program  may  refer  to  a  subroutine  subprogram  by  a 
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CiiLL  statement  which  specifies  the  name  of  the  subprogram  and  its 
arguments.  The  subroutine  subprogram  uses  one  or  more  of  its 
arguments  to  return  resialts.   Therefore,  the  arguments  that  are 
used  for  this  purpose  must  appear  on  the  left-hand  side  of  an 
arithmetic  statement  in  the  subprogram  or  in  an  input  list  within 
the  subprogram.         ■  '' 

The  correspondence  between  arguments  applys  in  the  case  of 
the  subroutine  subprogram  just  as  it  did  for  the  function  subpro- 
gram.  That  is,  the  arguments  listed  in  the  CALL  statement  must 
correspond  in  number,  order,  mode,  and  dimension. 

The  CALL  statement  refers  only  to  the  subroutine  subprogram 
whereas  the  BSTUEl-I  statement  is  used  by  both  the  function  and  sub- 
routine subprograms.   The  general  form  of  the  CALL  statement  is 
as  follows: 

CALL  "Name"  (al,  a2,  ....  ,  an)  . 

Name  is  the  name  of  the  subroutine  subprogram  being  called  and 

(al,  a2 ,  an)  are  the  arguments.   Each  argument  may  be  a 

fixed  or  floating  point  constant  or  variable  (vjlth  or  without  sub- 
scripts) ,  or  an  expression. 

A  priority  list  of  the  functions  in  order  of  their  power  and 
flexibility  would  be  as  follovfs: 

1.  Non-relocatable  Functions         ..  - 

2.  Library  Functions 

3.  Subroutine  Subprograms 
k.     Function  Subprograms 

5.  Arithmetic  Statement  Functions, 
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This  list  would  have  the  same  priority  of  difficulty  of  pro- 
gramming, the  first  being  difficult  and  the  last  being  norraal  For- 
tran programming. 

A  symbol  table  is  a  straightforward  list  of  all  encountered 
symbols  such  as  subroutine  names,  variable  names,  and  statement 
numbers  (Leeson,  1962).   A  symbol  table  look-up  operation  occurs 
each  time  a  constant,  variable,  or  a  statement  number  is  encountered. 
If  the  symbol  is  in  the  table  its  object  time  address  is  found  in 
the  corresponding  "Table  of  Addresses  of  Encountered  Symbols"  and 
is  used  in  the  generated  object  prograj^i.   If  the  syiabol  is  not  in 
the  table  it  is  placed  there  and  its  object  time  address  is  deter- 
mined and  stored  in  the  corresponding  address  table.   The  address 
is  then  used  as  above. 

A  brute  force  comparison  is  used  to  determine  if  a  symbol  is 
in  the  symbol  table  until  a  successful  comparison  is  made,  the 
temporary  end  of  the  symbol  table  is  found,  or  the  symbol  table 
is  found  to  be  full.   This  way,  the  mere  mention  of  a  symbol  de- 
fines it. 

The  first  symbols  in  the  symbol  table  are  the  names  and 
alternate  names  of  the  library  functions. 

With  the  use  of  the  symbol  table  and  the  table  of  encountered 
addresses,  Fortran  II  generates  a  series  of  branch  and  transmit 
instructions,  of  the  regular,  immediate,  and  floating  varieties, 
for  the  mainline  program.   The  P  address  depends  upon  the  operation 
desired  and  directs  control  to  a  subroutine.   The  Q,   address  depends 
upon  which  symbol  wants  action  at  that  time.  Both  will  come  from 
the  table  of  addresses  with  the  exception  of  operations  requiring 
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the  use  of  non- relocatable  subroutines.  These  have  a  fixed  lo- 
cation and  their  addresses  may  be  found  more  directly.   The  symbol 
table  and  the  table  of  addresses  are  both  compiling  aids  and  there- 
fore are  not  used  in  an  object  program. 

If  the  basic  processor  is  changed  we  lose  the  ability  of  in- 
corporating revisions  furnished  by  IBM,   If  we  work  at  the  library 
function  level  or  higher  with  very  little  additional  worl;;  our  devel- 
oped subroutines  may  be  used  by  a  different  processor  (on  a  machine 
other  than  the  1620).   For  this  reason,  alterations  and  additions 
irill  be  made  in  the  form  of  library  functions. 
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DECISION  SYSTEMS 

A  deoislon  system  is  the  algorithm  for  relating  interacting 
variables  in  a  problem  solution.  Deoislon  systems  can  be  simpli- 
fied and  put  in  a  losical  flow.   For  this  discussion,  then,  a 
decision  system  will  be  described  as  a  locical  and  orderly  method 
by  vihich  decisions  may  be  made. 

The  inefficiencies  of  a  decision  tree  have  already  been  dem- 
onstrated in  the  discussion  of  the  "IF"  statement.   For  simplifica- 
tion, this  discussion  will  be  reduced  to  one  of  "limited  entry",  or 
binary  decisions.  Even  when  using  "IP"  statements  in  Fortran,  the 
two-way  decision  is  often  adequate. 

The  decision  tree,  for  limited  entry  decisions,  is  sho^m  in 
figure  2,  vrhere  J   and  K  stand  for  Yes  and  no  respectively.  The 
use  of  two  decision  parameters  yields  four  paths  and  requires 
three  "IF"  statements  to  accurately  describe  the  system.   Thus 
20  parameters  yield  1,048,576  possible  paths  and  requires  1,048,575 
"IF"  statements  to  describe  the  system. 

An  example  of  a  decision  system  is  the  problem  of  credit  ap- 
proval (Kirk,  1965).  In  this  case  three  parameters  are  chosen  on 
which  to  base  a  decision.   They  are: 

1.  Is  the  credit  limit  OK? 

2.  Is  the  pay  experience  favorable? 

3.  Has  special  clearance  'oeen   obtained? 

The  decision  tree  to  describe  these  parameters  is  shovrn  in 
figure  2.   The  three  parameters  yield  eight  possible  branches  and 
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require  seven  comparisons  to  describe  the  systen. 

Structure  tables  are  an  advantageous  method  for  unambiguously 
describing  complex,  multi-variable,  multi-rule  decision  systems 
(Schmidt,  is6k) ,      Saoh  table  is  a  precise  statement  of  the  logical 
and  quantitative  relationships  supporting  a  particular  elementary 
process.  They  are  developed  in  terms  of  the  criteria  or  parameters 
affecting  the  problem  and  the  various  outcomes  which  may  result. 
An  example  of  a  decision  table  incorporating  the  above  parameters 
Is  described  in  figure  J, 

This  table  may  be  shortened  by  introducing  another  symbol 
other  than  the  Y  or  K,   This  symbol  is  "-"  and  means  "not  perti- 
nent",  Wienever  this  symbol  is  used,  the  parameter  in  question 
(the  row)  has  no  significance  in  that  rule  (the  oolu:.Tn) ,   Using 
this  new  symbol  the  decision  table  in  figure  3  is  reduced  to  figure 
i>.  Using  this  new  decision  table  we  E.ay  shorten  the  tree  in  figure 
2  to  that  shoTni  in  figure  5.   This  shortens  the  problem  to  one  of 
four  paths  and  three  "IP"  statements  by  eliminatins  redundancies. 

Another  economy  in  the  decision  table  is  gained  by  introduc- 
ing the  concept  of  "ELSE",   This  concept  merely  states:  if  no  rule 
in  a  decision  table  can  be  matched  then  the  action  specified  by 
"ELSE"  will  be  executed.   Using  this  idea  we  may  reduce  the  decision 
table  above  to  two  rules,  one  a  formal  rule  and  the  second  an  Implied 
rule  (figure  6). 

The  conversion  of  a  decision  table  into  a  decision  tree  is 
not  always  as  easy  as  above.   The  usual  case  has  been  to  construct 
the  general  tree  vrlth  all  of  its  redundancies  and  inefficiencies. 
Several  algorithms  by  which  a  fairly  consistant  and  efficient  tree 
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may  be  constructed  have  been  devised.   One  such  alsoritlni.  will  be 
presented  here,  This  algorithm  minlml2es  computer  storage  space 
required  for  the  resultant  progran.   It  is  a.lso  designed  to  pin- 
point any  contradictions  or  redundacies  among  the  rules  in  a  table 
(Pollaolc,  1965). 

■  Figure  ?  presents  an  example  of  a  limited-entry  decision  table 
that  could  easily  present  redundancies  and  inefficiencies  in  the 
normal  procedure  of  constructing  a  tree.   The  il-^   are  the  actions, 
the  Ci  are  the  conditions  or  parameters,  and  the  3]_  are  the  rules. 
Before  converting  the  decision  table  to  individual  comparisons  and 
the  series  of  branches  associated  with  each  path,  the  number  of 
Trritten  decision  rules  should  be  reduced  to  a  miniaiua  to  simplify 
the  use  of  the  algorithm.  As  an  example,  figure  7  may  be  reduced 
to  the  table  described  by  figure  8. 

It  would  be  helpful  to  present  a  general  description  of  hovi 
to  convert  the  decision  table  to  tree  form  before  the  algorithm  due 
to  Pollacl:  is  presented.  The  procedure  is  as  follows. 

One  row  of  the  original  decision  table  is  selected.   The 
criterion  for  selection  is  given  in  the  algorithm.   The  condition 
in  that  row  becomes  the  first  comparison  of  the  flowchart.   The 
original  decision  table  is  then  decomposed  into  two  subtables 
(containing  one  less  rov:) ,  or  a  subtable  and  a  rule,  or  only  two 
rules;  and  each  of  these  is  associated  with  each  branch  of  the 
comparison. 

This  is  continued  with  each  subtable  until  the  branches  lead 
only  to  completed  rules  or  "ELSE",  or  until  a  subtable  indicates 

that  the  original  table  contained  redundant  or  contradictory  rules. 

t-.  '■  ■     ^       ■     '   •    ■    I,- .  -       .:..■'■_.■' 
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The  alsorlthn  is  presented  without  proof,  hotrever  when  applied,  it 
has  proven  to  be  effective  (Pollaclc,  I965).   The  o-ojective  of  the 
algorithm  is  to  convert  a  decision  table  to  a  computer  program  and 
have  this  program  use  the  minimum  number  of  stora,",-e  locations,  and 
still  be  relatively  fast. 

Step  one  of  the  algorithin  is  a  check  for  redundancy  or  con- 
tradiction.  If  at  any  stage,  a  pair  of  rules  does  not  contain  at 
least  one  (Y,  S)  pair  in  any  of  its  roij,  redundancy  or  contradiction 
exists.   Such  is  the  case  for  rules  1  and  2  in  figure  9.   If  the 
actions  for  rules  1  and  2  are  identical  a  redundancy  exists  and 
rule  2  may  be  eliminated.   Figure  10  presents  an  example  of  this. 
Vlhere  the  actions  for  rule  1  differ  from  those  of  rules  2  and  3, 
a  contradiction  exists.   Where  the  actions  of  rule  1  are  identical 
with  those  of  rules  2  and  3,   a  redundancy  exists. 

Step  two  of  the  algorithm  is  to  malce  a  dash  count  and  deter- 
mine those  rows  which  have  a  minimum  dash  count,  A  dash  that  ap- 
pears in  a  rule  (column)  that  contains  r  dashes  is  counted  as  2 

r 
dashes.   For  each  rule,  2  is  denoted  as  the  column  count.   In 

each  row,  the  s'om  of  the  column  counts  corresponding  to  the  dashes 
in  that  row  is  called  the  dash  count.   A  row  dash  count  is  the  sum 
of  the  column  counts  of  those  jrules  that  have  dash  entries  in  the 
row.   This  is  illustrated  in  figure  11  in  which  row  2  has  the  mini- 
mum dash  count. 

Step  three  is  used  in  the  event  that  more  than  one  row  has  a 
minimum  dash  count.   It  is  to  select  that  row  which  has  a  maximum 
delta,  which  is  the  absolute  value  of  the  difference  of  the  X- 
oount  and  the  N-count.   The  Y-count  is  the  sum  of  the  column- 
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counts  correspondins  to  the  Y's  in  tho  roiT,   The  H-count  is  simi- 
lar.  Figure  12  is  an  example  of  a  decision  table  with  a  miniEum 
dash  count  in  all  of  the  rows.   Since  Cj^  had  the  maximum  delta, 
it  was  selected.   The  selected  row  is  called  k.      In  figure  12 
Ci=C,j., 

Step  four  is  to  discriminate  on  the  condition  in  row  !c.   This 
discrimination  has  two  branches,  each  of  which  leads  to  a  subtable 
X'Thlch  contains  one  or  more  rules,  with  one  row  less  than  the  orig- 
inal table  (row  k  is  deleted).   The  Y-branch  had  a  X  in  row  It. 
The  N-branoh  worlcs  similarly.   In  addition  both  subtables  contain 
those  rules  that  had  a  dash  in  row  k   of  the  previous  table. 

Step  five  is  to  go  back  to  step  one  if  the  subtable  of  inter- 
est contains  more  than  one  rule. 

Step  six  may  be  divided  into  four  smaller  steps  each  of  vrhioh 
handles  a  different  situation.   They  are  as  follows: 

a.  If  a  branch  leads  to  a  subtable  containing  one  rule, 
and  if  that  rule  contains  all  dashes,  then  replace 
the  subtable  with  the  rule  itself, 

b.  If  a  branch  leads  to  a  subtable  with  one  rule,  and  that 
rule  does  not  have  all  dashes  but  has  more  than  one  Y, 
K,  or  a  combination  of  Y  or  N,  choose  as  row  Ic  any  row 
that  has  no  dash.  The  selected  branch  will  indicate  a 
subtable  with  one  less  row  in  it.  The  opposing  branch 
will  Indicate  "ELSE", 

c.  If  a  branch  does  not  lead  to  a  subtable,  it  leads  to 
"ELSS" , 

d.  If  a  branch  leads  to  a  subtable  containing  only  one  rule. 
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and  that  rule  contains  only  one  oondltion  whose  value  is 
Y  or  K,  then  one  branch  of  the  discrimination  on  that 
condition  leads  to  the  rule,  the  other  branch  leads  to 
"ELSE". 
As  an  example  of  the  use  of  this  algoritlim,  refer  to  the 
table  in  figure  8.   The  application  of  the  algorithm  yields  figure 
13. 

Using  this  algorithm  one  may  reduce  a  decision  table  to  the 
minimum  tree  and  program  it  in  Fortran.   This  is  then  a  valuable 
programming  aid.   This  does  not,  however,  allow  a  programmer  to 
program  decision  tables  into  his  Fortran  program.   This  means,  also, 
that  the  programmer  must  recalculate  the  entire  decision  tree,  re- 
program  it,  and  reprocess  the  program  in  order  to  ma]:e  any  change 
in  the  decision  process,   A  possible  program  might  perform  the  pro- 
cess of  reducing  the  table  to  a  tree,  and  might  even  generate  the 
completed  "IF"  statements  on  cards.   The  compiling  process  woiad 
remain,  however.   It  would  be  more  desirable  to  change  the  decision 
table  at  object  level. 

It  is  possible  to  program  decision  tables  into  Fortran  with 
the  use  of  the  computed  "GO  TO"  statement  and  a  simple  arithmetic 
statement  (Veinott,  I966) ,   The  computed  "GO  TO"  statement  Indicates 
the  statement  that  is  to  be  executed  next.   However,  the  statement 
number  that  the  program  control  is  transferred  to  can  be  altered 
during  the  program.   The  general  form  of  the  computed  "GO  TO"  state- 
ment is 


GO  TO  (Kl,JJ2,...,Km),I 
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where  Nl,  K2,.,.,  Hm  are  statement  numbers  and  I  is  a  non-sub- 
scripted fixed  point  variable.   This  statement  causes  transfer  of 
control  to  the  first,  second,  third,  etc.,  statement  in  the  list 
depending  on  vfhether  the  value  of  I  is  1,  2,  3,  etc. 

The  approach  taken  considers  a  decision  table  as  a  multiple 
branch  within  a  program.   The  procedure  is  to  calculate  a  unique 
number  for  each  possible  set  of  conditions.   The  unique  numbers 
must  be  an  unbroken  series  of  consecutive  numbers  in  order  to  be 
used  as  a  branching  variable.   Figure  ^-   will  serve  as  an  example. 

A  new  table  is  constructed,  adding  a  "value"  column.   Since  there 

3 
are  2  =8  possible  combinations  of  three  conditions,  8  columns  will 

be  provided,  one  for  each  possible  combination.   The  8  columns  are 

numbered  from  0  to  7  inclusive.   X's  are  placed  in  the  ooliomns  so 

that  the  corresponding  "values"  add  to  the  number  heading  the 

colujLn.   This  is  represented  in  figure  14.   This  provides  a  unique 

number  for  each  possible  combination  with  a  consecutive  order. 

Since  the  series  contains  a  zero  it  is  necessary  to  add  one  to  make 

it  a  branching  variable.   Ordinarily  yes  is  denoted  by  1  and  no  Is 

denoted  by  0.   This  convention  is  carried  over  in  this  case.   In 

the  above  table  the  following  notation  is  used: 

11=1,  Credit  Limit  OK 

=0,  Otherwise 

12=1,  Pay  Experience  OK 

=0,  Otherwise 

13=1)  Special  Clearance 

=0,  Otherwise 
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Kl=  Statement  number  to  initiate  the  action 

"do  not  approve  order" 
N2=  Statenent  number  to  Initiate  the  action 

"approve  order" 

The  Fortran  program  for  this  table  is  as  follows: 

JUMP=l+Il+2*I2+ij-*I3 

GOTO  (N1,N2,N2,K2,N2,N2,N2,IJ2)  ,JUMP 

It  may  be  desirable  to  represent  one  or  more  conditions  by 
more  than  two  states.  Again,  let  there  be  m  conditions,  the  states 
of  each  of  which  is  indicated  by  the  value  of  variables  11,12,... 
,Im.   Let  the  various  conditions  have  Kl,K2,...,Km  mutually  exclu- 
sive states.   That  is,  the  conditions  themselves  are  represented  by 
the  "I"  variables:   each  of  these  "I"  variables  can  take  on  different 
values,  starting  with  zero,  to  express  the  state  of  this  particular 
condition.   The  number  of  states  of  any  condition. 

Since  the  states,  for  any  condition,  are  mutually  exclusive, 
by  definition,  only  one  state  can  exist  at  a  time  for  any  given 
condition.  The  number  of  combinations  or  "rules"  that  exist  will 
be: 

Number  of  rules  =  (Kl)  (K2) . . . (Km)=E.  (Veinott,  I966) 

For  convenience,  let  KNL   equal  the  number  of  states  of  the  nezt- 
to-the-last  condition. 

The  procedure  in  programming  such  a  table  is  to  set  up  H 
rules  and  identify  each  combination.   To  each  combination  a  state- 
ment number  must  be  assigned.   The  Fortran  program  would  be: 
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JOTiP=l+Il+ia*I2+Kl«K2*I3+ +(ia*K2, .  .*:fflFL)*Im 

GOTO    (N1,N2 ,Kr),JUKP 

For  an  example,  see  figure  15. 

This  technique  requires  the  exclusion  of  the  "-"  or  "not 
pertinent"  element  as  well  as  the  use  of  "ELSE".   This  means  that 
a  table  must  be  expanded  to  cover  all  possible  combinations  of 
conditions.   For  the  limited  entry  decision  table  2^  rules  are 
required  for  n  decision  parameters.   Finally,  changes  in  the  de- 
cision table  cannot  be  made  on  the  object  level.   The  user  must 
recompile  his  program  to  nalce  any  such  changes. 

A  more  efficient  and  faster  method  is  desirable,  as  well  as 
one  which  will  allow  changes  on  the  object  level.   An  algorithm 
which  should  best  fit  these  requirements  was  presented  by  Kirlc 
in  January,  I965.   It  is  as  follows. 

The  first  step  is  to  prepare  a  binary  image  of  the  condition 
portion  of  the  decision  table.   This  image,  called  the  "table 
vector"  and  consisting  of  "rule  vectors"  is  shoxmi  in  figure  16  for 
the  credit-approval  table  in  figare  3. 

The  table  matrix  and  the  data  vector  alone  do  not  have  suffic- 
ient information  for  solving  the  problem  because  of  the  use  of  a 
zero  for  nonpertinent  conditions  in  the  table  matrix.   A  masicing 
matrix  is  used  to  screen  out  nonpertinent  conditions  from  the  data 
vector  prior  to  scanning  the  table  matrix.   The  masking  matrix  is 
produced  by  replacing  a  Y   and  an  K  with  a  1,  and  "-"  with  0. 
Therefore,  a  logical  multiplication  will  insure  a  zero  in  the 
element  of  the  data  vector  corresponding  with  a  nonpertinent 
element  in  the  original  table.   The  masking  matrix  for  the  credit 
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approval  table  Is  as  in  figure  1?,  •-"■''' 

A  table  of  logical  multiplication,  shoxm  in  figure  18  is  used 
to  multiply  the  data  vector,  elenent  by  element,  by  the  appropriate 
masking  vector  prior  to  scanning.   The  result  is  then  compared  to 
the  corresponding  table  vector.   If  a  match  is  found  the  resulting 
action  is  executed;  if  not  the  next  masking  vector  and  table  vector 
are  tried  until  the  table  is  exhausted.   If  the  table  is  exhausted 
then  "ELSE"  is  executed. 

Since  the  IBM  1620  computer  is  not  a  binary  computer,  the  use 
of  characters  is  necessary.   The  logic  used  in  the  programming  of 
this  algorithm  is  presented  in  a  subsequent  section  of  this  thesis. 
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RAI''!DOri  KUi'iSSES 

A  Monte  Carlo  calculation  makes  systematic  use  of  random  num- 
bers. The  simulation  is  provided  with  intei-nally  generated  data 
that  has  the  same  characteristics  as  the  actual  data  (Boi'iEian  and 
Fetter,  196I) .  A  schematic  diagram  of  a  process  might  consist  of  a 
jagged,  irregular  figure  representing  the  collection  of  facts  as 
observed.   An  idealized  model,  represented  by  a  smooth  curve  may  be 
derived  from  the  schematic  diagram  of  observations.   This  model, 
used  with  a  random  number  generator,  produces  an  artificially  irreg- 
ular figure,  representing  the  simulated  experience.   The  simulated 
data  is  used  in  mathematical  models  for  further  simulations.   In  ad- 
dition to  providing  more  extensive  data,  you  can  generate  your  data 
according  to  any  icnoifin  rule  or  give  it  any  special  characteristics 
that  you  wish  in  order  to  study  the  effects  on  your  simulation  model. 

The  scheme  choosen  for  generating  random  numbers  has  the  fol- 
loTfing  features:  (a)  the  numbers  generated  are  uniforialy  distributed 
betvreen  0  and  1;  (b)  there  is  no  serial  correlation  between  succes- 
sive numbers  in  the  sequence;  (c)  about  50  million  numbers  can  be 
generated  before  the  sequence  repeats  itself;  (d)  the  calculation  is 
easy  to  perform  on  an  electronic  computer  (Hull,  I962) . 

The  scheme  is  this:  Start  with  any  ten-digit  number  of  the  form 
xyzOOOOOOl;  call  it  Kg.  Thereafter 

Hjj  =  K  •  Rj^_i  (mod  10'''°) 

vrhere  IC  is  a  fixed  multiplier,    which  should  be  a  ten  digit   odd  power 
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of  a  prine  that  is  relatively  prime  to  10.   Possible  values  for  K 
include  (Hull,  1962), 

3^5  =  1,162,261,46? 
7II  =  1,977,326,743 
11^  =  2,357,947,691. 

The  recursion  relation  expressed  above  in  symbols  can  be  trans- 
lated: To  get  the  nest  nujuber  in  the  sequence,  multiply  the  previous 
number  by  K.  The  result  will  be  a  20  digit  number.   Take  the  last 
(right-hand)  ten  digits  of  the  product  as  the  next  number  in  the  se- 
quence.  Treat  the  number  as  a  ten-digit  decimal. 

This  scheme  will  generate  rectangular  or  uniform  random  numbers 
which  in  turn  may  be  used  to  ssxiple  from  any  distribution  desired. 
It  was  coded  in  the  Symbolic  Programming  System  (SPS)  language  for 
the  1620  and  compiled  as  a  library  function  for  Fortran  II  (see  page 

).   ';:.-,.:  ■       : 

Non-uniform  random  numbers  are  occasionally  needed  in  order  to 
more  properly  describe  a  particular  process.   For  one  dimensional 
distributions  we  need  only  to  solve  the  equation  x=?(y)  for  y,  vrhere 
X  is  uniformly  distributed,  and  where  P(y)  is  the  required  (cumula- 
tive) distribution  function.   For  example,  if  y  is  to  be  exponentially 
distributed,  the  following  distribution  function  is  used: 

X  =  1  -  e'^  and 
y  =  ln(l  -  x). 

An  arithmetic  statement  function  to  generate  random  n-jmbers  repre- 
senting an  exponential  distribution  could  be  as  follows. 
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E3CPB1I(X)  =  L0C-?(1.  -  EI.im(X)) 

The  use  of  these  functions  has  been  previously  described. 

Many  simulations  require  the  use  of  arbitrary  distributions, 

iThich  may  be  used  with  a  table  loo!c-up  operation.   For  a  detailed 

study  on  the  vise   of  arbitrary  distributions,  see  Starr  and  i.iller 

(12), 

Since  its  developenent,  the  random  number  generator  has  been 
tested  by  VJichlan  (l'^-).   ^iie  numbers  were  found  to  be  rcots.n^-ular 

on  the  five  per-oent  level.   The  Chi-squared  and  Serial  Correlation 

tests  were  used. 
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IKPLEJIENTATION  OF  TH3  KODIFISD  rHOCSSSOH 

The  Decision  Table  sutooutliie  h;r.s  been  coded  into  a  relocatable 
library  subroutine.   There  is  present  a  need  for  the  capability  to 
set  up  data  vectors,  i.e.  present  the  conditions  as  they  actually 
er.ist  in  the  fom  of  a  vector.   This  capability  must  be  in  the  form 
of  relocatable  subroutines  to  be  conpatible  with  the  Decision  Table 
subroutine.   A  one  parameter  data  vector,  with  the  help  of  the  proper 
subroutine,  could  be  set  up  -rith  the  followinc  Fortran  statement. 

DATA  =  EQ(T-S) 

If  T  and  S  are  equal,  a  1  will  becoae  the  entire  data  vector;  if  not 
a  0  will  be  produced.   The  functions  available  for  oonstractins  the 

data  vector  a.re; 

-  ■■  . .   ^_ 

SQ   Equal 

UK   Not  equal      '  "'■' 

GE   Greater  than 

LR   Less  than 

GE   Greater  than  or  equal         ,  ^  ^^^__' 

LE   Less  than  or  equal 

KC   No  change  (for  data  vectors  that  have  already  been 

produced) 
YES  Always  generates  a  1    ■ 
HO   Always  generates  a  0 
EDATA  The  entire  data  vector  will  be  read  in. 

An  example  of  the  use  of  these  functions  is  illustrated  by  the 
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following  Fortran  statement. 

DATA=EQ  { T-S )  +GH  ( U-T+LRU-S )  -i-GE  ( T-U)  +LE  ( U-T ) 
if  T=l,  S=2,  and  U=3  the  following  data  vector  will  be  produced. 

011000 

The  data  vector  may  then  be  tested  after  a  table  has  been  read 
into  core.   This  is  accomplished  by  the  following  statement. 

ACTI0N=BDTB(1) 

This  will  cause  a  decision  table  to  be  read  into  the  proper  place 
in  core  (common).   It  will  also  be  named  table  niimber  1.  There- 
after a  table  may  be  referred  to  by  its  number  in  testing  as 
follows: 

JUMP=TEST(1). 

The  last  data  vector  produced  will  then  be  tested  by  decision 
table  number  1.   The  resulting  value  will  be  assigned  to  the  branch- 
ing variable,  JUMP.   JUMP  is  then  used  in  a  computed  GO  TO  state- 
ment. 

The  credit  approval  table  will  serve  as  an  example  (figure  4), 
The  following  variables  will  be  used  in  the  example: 

JUMP  =  the  branching  variable 

OEDLMT  =  the  customers  credit  limit 

CUSDBT  =  the  amount  owed  by  the  customer 

OHDAMT  =  the  dollar  amount  of  the  order 

X  =  a  dummy  variable 
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The  Fortran  program  is  as   follows; 

DATA=GE(CRDLMT-(CUSDBT+OBDA]yiT)  )+YES(X)+NO(X) 

ACTI0N=HDTB(1) 

JUMP=TEST(1) 

GO  TO   (1,1, 1,2), JUMP 

The  DATA  statement  may  vary  and  options  may  be  placed  under 
switch  control.   When  DATA  is  executed  a  data  vector  is  produced 
and  stored  in  core  according  to  the  function  names  and  the  values 
of  their  arguments.   When  the  function  EDTB(l)  is  executed  a 
table  will  be  read  into  common  and  labeled  with  a  1.   In  this 
case  the  table  will  look  like  the  folloifing: 

Oif-03 
1— + 
01-* 
001* 

ooot 

vjhere  0403  means  that  there  are  four  rules  and  three  conditions 
in  the  table.   For  each  table  the  limit  on  conditions  is  79  and 
the  limit  on  rules  is  99.   The  rules  may  be  presented  in  any  order 
but  the  desired  actions  must  be  expressed  in  the  same  order  in  the 
computed  GO  TO  statement  used  in  the  source  program.   ACTION  and 
DATA  are  dummy  variables  but  must  be  floating  point  variables. 
The  data  vector  must  fit  the  table  it  is  intended  for.   The  decis- 
ion table  and  data  vector  must  be  based  on  the  same  number  of 
conditions  or  parameters.   In  addition,  the  data  vector  must  be 
used  before  another  data  vector  is  produced.   However,  a  decision 
system  may  contain  more  than  one  decision  table.  Tables  may  also 
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vary  in  length. 

When  the  arithmetic  statement  defining  JUMP  is  encountered 
a  value  will  be  assigned  to  the  variable,  JUMP.  Assuming  this 
particular  customer  has  exceeded  his  credit  limit,  the  value  of 
the  data  vector  is  as  follows. 

010 

The  resulting  value  of  JUMP  will  be  2  and  statement  number  1 
will  be  the  next  statement  executed  within  the  program. 

If  a  data  vector  is  to  be  read  in,  the  entire  data  vector 
must  be  ready  and  must  be  right  justified  on  the  data  card  with 
a  flag  over  the  left-most  element  on  the  card  as  follows. 

columns  78-80         010 

The  data  vector  should  have  the  same  number  of  conditions  as  does 
the  table  used  for  testing.    -r 

The  purpose  of  the  "data  arranging"  subroutines  is  to  set 
up  the  data  vectors  in  a  place  in  core  which  is  accessable  to  all 
of  the  relocatable  subroutines.   In  this  ■way   any  subroutine  can 
make  changes  to  or  use  the  data  vector  as  it  must.  The  purpose  of 
the  "table  reading  and  testing"  subroutine  is  to  read  a  table  into 
core,  label  it  and  set  up  the  table  and  masking  matrices.  At 
the  time  of  testing  the  logical  multiplication  is  performed  on 
the  data  vector  to  mask  it  and  the  result  is  compared  to  the  table 
vector.   Each  time  a  comparison  is  made  a  counter  is  incremented 
by  one.   I'Jhen  a  comparison  is  successfully  made  this  counter  be- 
comes the  value  of  the  branching  variable.   The  subroutines  and 
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their  flow-charts  are  presented  In  the  appendix. 

In  order  to  implement  relocatable  library  functions,  they 
must  first  be  coded  into  a  1620  SPS  source  prosram.  The  origin 
is  assigned  at  location  10,000.  All  P  and  Q  addresses  that  are 
relative  to  the  origin  must  be  labeled  as  such  by  placing  a  flag 
in  the  0-,    and  Op  positions  of  that  same  instruction,  respectively. 

The  resi;ilting  condensed  object,  w-ith  the  first  2  and  last  7 
cards  removed,  must  be  preceded  by  a  header  card  punched  vjith  the 
following  information: 

Columns  1-  5  XX30QC  Total  number  of  storage  locations 
required  by  the  subroutine.  This 
number  must  be  even. 

Columns  11-12  YY      The  alternate  subroutine  number 
if  any. 

Colianns  16-20  WVfWliW   The  alternate  entry  point  if  any. 

Column  63       -A  record  mark. 
_.  Columns  75-80  XXOOOl   Card  sequence  number,  where  XX 
is  the  subroutine  number. 

The  object  deck,  minus  the  header  card,  must  be  renumbered  start- 
ing with  XX0002,  where  XX  is  the  subroutine  number. 

If  a  subroutine  contains  two  entry  points  the  card  sequence 
numbering  does  not  change  and  is  done  as  if  no  alternate  entry 
vrere  in  effect.   The  information  given  on  the  header  is  suffi- 
cient for  the  use  of  the  alternate  entry  point.   The  packet  of 
cards  can  be  inserted  between  any  two  subroutines  in  the  library 
subroutine  deck. 
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It  is  also  necessary  to  increase  the  number  in  columns  1-2 
of  card  030OO  of  the  Pass  I  deck  to  reflect  a  new  total  number  of 
relocatable  subroutines.  Then  add  the  subroutine  name(s)  at  the 
end  of  the  last  card.   If  a  new  card  is  required,  it  must  be  num- 
bered consecutively. 

Linkage  into  the  subroutine  is  provided  with  one  of  the 
"Branch  and  Transmit"  instructions.   For  instance: 

BTM   SUBH,   A  .  >  '.  1 

will  branch  to  SUBE  and  carry  A  I'lith  it.   It  will  place  the  field 
represented  by  A  into  core  just  before  the  entry  point  of  SUBE. 
When  the  subroutine  is  finished  a  "Branch  Back"  will  be  encountered 
and  control  will  be  transferred  back  to  the  mainline  program. 
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TRIAL  PEOGEAJIS 

In  order  to  test  the  effectiveness  of  the  decision  processes 
and  Illustrate  the  use  of  the  random  number  generator,  a  system 
was  chosen  for  simtilation.   This  system  entails  a  pump  and  a  tank. 
The  pump  is  used  to  fill  the  tank  and  may  have  various  capacities. 
The  capacity  of  the  pump  is  stated  in  terms  of  the  depth  of  liquid 
replaced  in  the  tanlc  per  pass  through  the  simulation.   Thus,  a 
continuous  process  is  simplified  into  a  discrete  process.   The 
pvunp  will  also  have  a  tendency  to  heat  during  the  process  of  pump- 
ing.  It  is  assumed  that  the  process  of  pumping  will  increase  the 
temperature  of  the  pump  by  one  degree  during  one  pass  through  the 
simulation.   Initially  the  pump  is  70  degrees  and  the  upper  safety 
limit  is  180  degrees.   It  Is  assumed  that  the  atmosphere  is  70 
degrees. 

The  loss  of  heat  from  the  pump  depends  upon  the  temperature 
of  the  pump.  The  computer  program  removes  heat  from  the  pump  by 
finding  the  temperature  gradient  between  the  pump  and  the  atmos- 
phere and  subtracting  one  hundreth  of  it  from  the  pump  tempera- 
ture during  each  pass  through  the  simulation  when  the  pump  is 
inactive.  *_,',,.,  '..  ,  /• 

The  tank  has  a  valve  which  may  be  used  to  remove  liquid  from 
the  tank  at  the  rate  of  "B"  feet  per  pass  through  the  simulation. 
The  control  of  this  valve  Is  placed  linder  a  user  of  users.   The 
program  places  the  control  of  this  valve  under  a  switch  on  the 
1620  console  or  will  simulate  the  use  of  the  valve  with  the  ran- 
dom number  generator  at  the  discretion  of  the  operator.   It  is 
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considered  desirable  to  keep  the  level  of  the  tanli:  'between  four 
and  five  feet. 

In  addition,  two  alarms  are  provided.  The  first  is  a  heat 
alarm  and  will  be  turned  on  when  the  temperature  of  the  pump  ex- 
ceeds 180  degrees.   The  second  is  a  "tank  level  high"  alarm  and 
will  be  turned  on  when  the  liquid  level  in  the  tank  exceeds  five 
feet.   The  decision  table  used  to  govern  the  pump  and  the  alaims 
is  depicted  by  figure  I9. 

This  entire  process  has  been  simulated  in  four  different 
ways  each  of  which  differs  In  the  decision  process.   The  general 
flow  diagram  is  as  presented  in  figure  20.   They  are  as  follows: 

1.  Using  a  tree  constructed  by  eliminating  the  table  (figure 
19)  one  condition  at  a  time  In  the  order  of  their  occur- 
rence. 

2.  Using  a  tree  constructed  by  eliminating  the  table  accor- 
ding to  the  "dash  count"  of  the  conditions  (thus  construct- 
ing an  optimal  tree). 

3.  Programming  the  table  in  the  Fortran  language. 

1+,      Using  the  table  on  the  object  level,  thus  using  the  afore- 
mentioned library  subroutines. 

The  tree  which  resulted  from  reducing  the  table  one  condition 
at  a  time  is  illustrated  in  figure  21.   This  reduction  was  accom- 
plished by  first  discriminating  on  condition  number  one  (tank  in 
use)  and  progressing  directly  through  the  rest  of  the  conditions 
until  the  conversion  was  complete.   The  resulting  program  is 
Illustrated  by  figure  22. 

The  tree  which  resulted  from  the  reduction  of  the  table 
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according  to  the  dash  count  is  as  described  by  figure  23.   Accord- 
ing to  the  proponents  of  the  algorithm  used,  this  is  the  smallest 
tree  that  can  be  produced.   The  resulting  program  is  illustrated 
by  figure  24. 

Programming  the  decision  table  in  the  Fortran  language  is 
accomplished  by  expanding  the  decision  table  into  the  form  indi- 
cated in  figure  25.  This  is  then  programmed  with  the  use  of  a 
single  arithmetic  statement  and  a  computed  "GO  TO"  statement. 
The  resulting  program  is  illustrated  by  figure  26. 

The  use  of  the  decision  table  on  the  object  level  calls 
for  a  reduction  for  the  decision  table  as  depicted  in  figure  2?. 
The  resulting  program  is  illustrated  by  fiugre  28. 

The  resulting  programs  are  compared  in  figure  29.  This  com- 
parison is  made  to  illustrate  the  difference  in  programming,  time 
and  core  storage  necessary  to  perform  the  simulation  with  1000 
passes.  Although  the  tank  use  was  controlled  by  the  random  number 
generator,  the  same  tank  usage  was  experienced  since  the  same  ar- 
gument controlled  the  random  number  gnerator  in  all  of  the  programs 
in  the  comparison. 

The  fastest  programs  were  the  programs  which  used  the  decision 
trees.  Both  had  times  of  1.901  minutes.   The  slowest  program  was 
the  program  using  the  SPS  subroutines  with  a  time  of  5.222  minutes. 
The  shortest  mainline  program  is  the  program  using  the  SPS  sub- 
routines.  However,  when  the  subroutines  are  taken  into  considera- 
tion this  program  becomes  the  longest.   The  shortest  program  (in- 
cluding the  required  subroutines)  is  the  decision  tree  produced  by 
Pollack's  algorithm.  However,  the  subroutines  are  of  a  constant 
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length  and  the  mainline  programs  are  of  a  variable  length.   For 
programs  requiring  more  complicated  decision  systems  the  program 
produced  with  Pollack's  algorithm  would  soon  prove  to  be  shorter. 
There  would  be  a  break-even  point  where  the  length  of  the  improved 
decision  tree  program  would  be  the  same  as  the  length  of  a  program 
utilizing  the  decision  table  subroutines.   Any  larger  systems  would 
yield  a  shorter  program  with  these  subroutines. 


C  '"■ 
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CONCLUSION 

The  comparison  of  the  different  trial  programs  yields  the 
following  information. 

1)  Small  to  medium  decision  systems  yield  shorter  programs 

2 

when  reduced  to  a  tree  with  Pollack's  algorithm  . 

2)  Larger  decision  systems  would  probably  yield  the  short- 
est programs  using  decision  tables  and  the  special  sub- 
routines. 

3)  The  particular  example  shoirn  above  indicates  that  faster 
programs  are  feasible,  using  a  decision  tree.  There  is, 
however,  nothing  on  which  to  base  a  generalization  about 
the  speed  of  resulting  programs. 

In  addition,  the  SPS  subroutines  allow  changes  to  be  made  on  the 
object  level  of  compiled  programs,  thus  allowing  a  greater  flexi- 
bility. This  is  important  for  simulations  because  you  can  try 
several  different  rules  and  arrangements  and  study  the  effects 
they  have  on  an  organization  without  recompiling  any  programs. 
The  SPS  subroutines  are  of  the  following  lengths. 

Eandom  Number  Generator 240 

Bead  Table  and  Test 2164 

Greater  than  or  equal  and  Less  than  equal 290 

Greater  than  and  Less  than 290 

Equal  and  Unequal 280 


^Pollack  has  met  with  recent  criticism  of  his  algorithm.   His 
assumption  has  been  proven  not  infallable  by  Sprague  (11). 
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Xes  and.  No 19^ 

No  Change 90 

Read  Data 36 

Some  of  these  are  quite  small  and  efi-ioient;  others  (due  to  their 
complexity)  are  large  and  bulky.   It  is  quite  probable  that  some 
of  the  subroutines  oould  be  shortened  and/or  made  to  operate  fast- 
er by  another  approach,  either  to  the  problem  solution  or  to  the 
techniques  used  in  programming  the  solution.   Due  to  its  length 
and  bulk  the  "Head  Table  and  Test"  subroutine  could  probably  be 
shortened  considerably. 


FIGURES 
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PLATE  I 
Figure  1  is  an  exsxmle   of  a  decision  tree  in  which  each 
condition  has  three  states.   Figure  2  illustrates  a  binary 
or  limited-entry  decision  tree. 
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PLATS  II 
Figure  3  is  a  credit  approval  decision  talDle  which  is 
expanded  to  include  all  possible  oases.   Figure  k   is  the  ssune 
table  after  the  rules  have  been  combined  with  "not-pertinent" 
elements.   Figure  5  is  the  credit  approval  table  illustrated 
in  the  form  of  a  tree. 
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-   ■        ,--,   PLATS  III 
Figure  6  is  the  credit  approval  table  simplif ied  by 
using  the  implied  decision  rule  -  "ELSE".  Fisure  7  is  any 
decision  table.  Pigure  8  is  the  simplified  version  of  figure 
7  after  the  introduction  of  the  "not-pertinent"  element. 
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PLATE  IV 
Figures  9  and  10  illtistrate  decision  tables  which  contain 
redundant  or  Inoonsistant  rules.  Fisure  11  is  a  decision 
table  which  Indicates  the  use  of  a  dash  count.   Figure  12 
illustrates  the  use  of  a  dash  count  and  the  delta  of  each 
condition. 
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PLATE  V 
figure  13  Illustrates  the  application  of  Pollaolc's  al- 
gorithm for  the  expansion  of  decision  tables  into  mininum 
decision  trees. 
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PLAT3  VI  ■-  '  *-   -^   ■  '  - 
Figure  Ik   illustrates  the  credit  approval  table  modified 
for  prograimiiinc  In  Fortran.  Pl^are  15  Illustrates  a  table 
containing  conditions  with  more  than  two  states  modified  for 

Fortran  procrsjiiniing. 
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PLATE  XII   ,■-,. 
Figure  16  is  the  table  matrix  for  the  credit  approval 
table.   Figure  1?  is  the  maslcing  matrix  i''or  the  credit  ap- 
proval table.   Figure  18  illustrates  a  table  of  logical  mul- 
tiplication. 
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PLATS  VTTI 
Fisiire  19  is  the  decision  table  for  the  "pump  s.nd  tan!:" 
simulation. 
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Figure  20  is  the  flow  logic  used  to  simulate  the  pump         ,j 
and  tanli:.  ; 
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Pigiire  21  is  the  uninproved  decision  tree  for  the  p-LEip 
simulation  taole. 
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PLATS  XI  ■  ■■  " 

?lsure  22  (3  pa^es)  is  the  program  for  the  unimproved 
decision  tree  for  the  punp  simulation. 
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■"iGure  22   continued. 
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Figure  22   concluded. 
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PLAT3  XIV 
Plsure  23  is  the  Improved  decision  tree  for  the  pimp 
simulation. 
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PLATS  ;CV 
Figure  Z^■   is  the  Fortran  program  for  the  Improved  decis- 
ion tree  for  the  pump  simulation. 
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Figure  25  is  the  expanded  decision  table  for  prosram- 
ming  the  puap  simulation  decision  table  in  Fortran, 
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Figure   26   continued. 
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Figure  28  Is  the  prosraa  for  simula-tln^  the  puiap  and. 
tanlc  with  the  Kodifled  processor. 
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Figure  28  continued. 
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Figure  28  concluded.   Decision  table  input  for  the  modified 
processor.  Figure  29  is  a  comparison  of  the  four  different  pump 
simulations. 
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240 

14112 

Program  U- 

5,222 

13514 

3066 

I658O 

Fig.  29 
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APPENDIX  1 
Suliroutines 

1)  Random  Niaaber  Generator 

Figure  Al Flowchart 

Program  Al SPS  Source  Program 

2)  Read  Table  and  Decision  Subroutine 

Figure  A2(a) , Head  Table 

Figure  A2Cb) Set  Up  Matrices 

Figure  A3  (a) Determine  Branching  Variable 

Figure  A3(b) Mask  Data  Vector 

Program  A2 Read  Table  and  Decision  Subroutine 

3)  Set  Up  Data  Vector  Subroutines 

Program  A3... Greater  Than  or  Equal,  Less  Than  or  Equal 

Program  A^J- Greater  Than,  Less  Than 

Program  A5 Equal,  Unequal 

Program  a6 Yes,  No 

Program  A? No  Change 

Program  A8 Read  Data 
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The  objective  of  this  thesis  was  to  develope  a  system  for 
the  IBM  1620  computer  which  would  facilitate  the  simulation  of 
systems.   The  Fortran  II  programming  language  was  chosen  as  a 
basis  for  further  study  and  implementation.  Four  methods  for 
programming  decision  systems  were  studied  and  implemented  in 
Fortran  either  as  Fortran  programs  or  as  subroutines  to  Fortran. 
The  four  methods  were  each  used  to  program  a  simulation  and  the 
results  were  compared.  The  four  methods  xrere;   normal  Fortran 
Programming  (using  a  decision  tree),  an  algorithm  for  minimum 
Fortran  programming  (improved  decision  tree),  a  decision  tables 
programmed  in  Fortran,  and  decision  tables  used  with  Fortran  ob- 
ject programs  (using  SPS  subroutines). 

The  results  indicate  that  shorter  and  probably  faster  pro- 
grams can  be  -written  using  the  minimum  decision  tree  programmed 
in  Fortran  for  small  to  medium  decision  systems.   However,  for 
large  decision  systems,  shorter  and  possible  faster  programs  may 
be  vjritten  if  use  is  made  of  the  SPS  subroutines. 


